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ABSTRACT 

The Extragalactic Lens VLBI Imaging Survey (ELVIS) searches for central images of lensed radio 
quasars, in order to measure the central density profiles of distant galaxies. Here we present sensitive 
multi-epoch Very Long Baseline Array (VLBA) observations of PMN J1838-3427 at 8 GHz, with a 1 a 
noise level of 38 /LtJy beam -1 . Based on the absence of a central image of the background source at this 
level, we explore the possibilities for the central matter distribution in the lens galaxy. A power-law 
density profile, p cx r~ 7 , must have 7 > 1.93. Thus the density profile is close to an isothermal profile 
(7 = 2) or steeper. The upper limit on any constant-density core in an otherwise isothermal profile 
is < 5 parsecs. We also derive the constraints on models in which the density profile is isothermal on 
kiloparsec scales, but is allowed to have a different power law in the central ~ 100 parsecs. If the lens 
galaxy harbors a supermassive black hole, the galaxy profile is allowed to be shallower, but for the 
expected black hole mass the galaxy profile must still be close to isothermal or steeper. 
Subject headings: gravitational lensing — galaxies: elliptical and lenticular, cD 



1. INTRODUCTION 

The central regions of galaxies (r < 1 kpc) are a topic 
of great interest, but one for which it is difficult to com- 
pare the results of observations and theoretical simula- 
tions. The observational difficulty is especially severe 
for galaxies at cosmological distances. The Extragalactic 
Lens VLBI Imaging Survey (ELVIS) is motivated mainly 
by the desire to measure (or place interesting constraints 
on) the central matter distribution in galaxies at signifi- 
cant redshift, by searching for central images of gravita- 
tional lenses. As this is the first ELVIS publication, we 
elaborate on the theoretical and observational situation. 

On the theoretical side, cosmological dark matter sim- 
ulations produce dark matter halos with a universal den- 
sity profile that goes as p oc r~ 3 at large radii and 
p oc r _7 ,7 = 1 — 1.5 at small radii (Navarro et al. 
1997; Moore et al. 1999). The transition occurs at a 
scale radius r s , typically 10 — 30 kpc for halos of galac- 
tic mass. The main difficulty with interpreting the re- 
sults for galaxy interiors is that on scales ~ 10 kpc and 
smaller, the baryonic component is expected to modify 
the dark matter halo significantly. Adiabatic contraction 
(Blumcnthal et al. 1986) and similar models (Sellwood 
& McGaugh 2005) model the contraction of dark matter 
under the gravitation of a baryonic disk or bulge. The to- 
tal matter profile steepens and becomes close to isother- 
mal (p oc r™ 2 ) on scales of a few kpc (Kazantzidis et al. 
2006). Similar results are found from hydrodynamic sim- 
ulations that use cooling, gas dynamics, star formation 
and gravitation to model the dark matter and baryons 
within individual halos. The baryonic matter dissipates 
energy and collects at the halo center, becoming the dom- 
inant component inside radii 1.5 — 5 kpc (Gnedin et al. 
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2004; Maccio et al. 2006). The dark matter contracts 
inwards under the gravitational influence of the baryons, 
again giving a total matter profile that is close to isother- 
mal. At present the smallest scales probed by the hydro- 
dynamic simulations are 0.3-1 kpc. 

On the observational side, the density distributions of 
massive galaxies can be directly probed through dynami- 
cal studies, at least in the nearby universe. Gerhard et al. 
(2001) have modelled the dynamics of a sample of large 
early- type galaxies, using photometry and kinematic line 
profiles. The dark matter fraction is only 10—40% within 
the effective radius i? e , while the rotation curves are flat 
on scales larger than 0.2i? e , indicating an isothermal den- 
sity profile on these scales. Typically R e — 4 — 10 kpc 
(Kronawitter et al. 2000), so these results agree with the 
simulations described above. The density profile of an 
early type galaxy is isothermal at radii of a few kpc, and 
baryons represent most of the mass inside this radius. 

Projected surface brightness profiles of nearby galaxies 
(z = 0.002 — 0.005) have been observed with the Hubble 
Space Telescope (Lauer, T. R. et al. 1995). The angular 
resolution is O'.'l, which corresponds to a physical size 
of ~ 10 pc for the typical galaxy in this sample. The 
surface brightness profiles are well fit by a broken power 
law, with steep outer exponents, shallower inner expo- 
nents, and break radii between 10 and 1000 pc. Based 
on their inner profiles I(R) oc i? _/3 the galaxies can be 
classified into two populations: steep cusps with 8 ~ 1 
and flatter cores with 8 = — 0.3. A surface density 
power law E(i£) oc R~P corresponds to a density power 
law p(r) oc r -7 , where 7 = 8 + 1. Thus the luminous 
matter in a cuspy galaxy has an isothermal distribution 
(£(i£) oc i? _1 ) to within ~ 10 pc of the galaxy's center, 
while the luminous matter in a core galaxy breaks to a 
shallower profile at some radius < 1 kpc. This is a good 
approximation to the total mass profile, as the stars seem 
to represent most of the mass at these radii. 

Gravitational lenses provide information on the mass 
profiles of more distant (and therefore younger) galaxies. 
For a lens galaxy at z = 0.3 — 1.0, the relative positions 
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and fluxes of the bright lensed images of a background 
source constrain the matter profile interior to a few kpc 
from the lens galaxy center. Detailed studies of <~ 20 
gravitational lenses find that the early type lens galaxies 
have density profiles which are very close to isothermal 
on these scales (Rusin et al. 2003; Koopmans et al. 2006). 
Distant early-type galaxies have similar profiles to those 
nearby: they are isothermal at galactic radii of a few kpc. 

But what about the central few hundred parsecs? 
Here, too, gravitational lenses can help, through the 
properties of the "central image." In theory, a non- 
singular galactic profile produces an odd number of im- 
ages (Dyer & Roeder 1980; Burke 1981). One image 
forms near the center of the lens galaxy, where it is ex- 
pected to be highly demagnified by the large surface den- 
sity at that position (Narasimha et al. 1986). Due to 
the demagnification, the faint central image is rarely ob- 
served, leaving two or four bright images. In cases where 
the density profile is singular, with a central cusp that 
is stronger than isothermal (p <~ r~ 7 with 7 > 2), no 
central image is produced even in theory. 

While they are hard to observe, central images probe 
the inner 10-100 pc of very distant galaxies. Winn et al. 
(2003, 2004) have confirmed the existence of a central 
image produced by an isolated lens galaxy in the lens 
PMN J1632-0033, using radio observations. We discuss 
this object and its matter profile further in Section 5. 
Inada, N. et al. (2005) found a central image generated 
by the combined profile of a cluster and a massive galaxy, 
using the Sloan Digital Sky Survey and the Hubble Space 
Telescope (HST) . The mass modelling here is dependent 
on both the galaxy and the cluster profile. 

Most galaxies host a super-massive black hole (SMBH) 
whose mass correlates closely with the properties of 
the central stellar bulge (Kormendy & Richstone 1995; 
Magorrian, J. et al. 1998), particularly the velocity dis- 
persion (the M — a relation; Ferrarese & Merritt (2000): 
Gebhardt, K. et al. (2000); Tremaine, S. et al. (2002)). 
A SMBH in the lens galaxy affects the central image: it 
can destroy the central image, or split the central image 
into two images, one of which is directly attributable to 
the black hole (Mao et al. 2001; Bowman et al. 2004). In 
the latter case, the properties of the central-image pair 
could allow for the measurement of the black hole mass 
in an ordinary galaxy at significant redshift (Rusin et al. 
2005). 

Based on the surface-brightness profiles of nearby 
early-type galaxies measured with HST, Keeton (2003) 
predicted the distribution of core image magnifications. 
He found a broad distribution, from 10~ 45 to 10 _1 , with 
a most probable magnification of approximately 10~ 25 . 
Adding a SMBH had little effect on central-image de- 
tcctability: the magnification of the central image was 
strongly affected only when the magnification of the 
galaxy alone was already very low. Rusin & Ma (2001) 
placed limits on central mass distributions using radio 
observations of six lenses from the Cosmic Lens All-Sky 
Survey (CLASS; Browne, I. W. A. et al. (2003); Myers, 
S. T. et al. (2003)). Their basic result was that if the 
density profiles in the lens galaxies are taken to be power 
laws in radius, then they must be nearly isothermal or 
steeper. The constraints were somewhat weaker when a 
SMBH was included. 

This paper is organized as follows. Section 2 describes 



the design of our survey. Section 3 presents new obser- 
vations of our first target, PMN J1838-3427. Section 4 
presents models for the density profile in the lens galaxy 
of this system, and the final section summarizes the re- 
sults. 

2. ELVIS 

The Extragalactic Lens VLBI Imaging Survey (ELVIS) 
will involve sensitive, high-angular-resolution radio ob- 
servations of many of the known cases of gravitational 
lensing of a radio-loud quasar. The traditional reasons to 
conduct such observations are to confirm cases of gravita- 
tional lensing, and to observe correspondences between 
lensed radio jets in order to refine models of the lens 
galaxy. ELVIS is the first survey (to our knowledge) mo- 
tivated by the search for central images. As such, our 
highest-priority targets are those that are most favor- 
able for central-image hunting: radio-loud, asymmetric 
two-image lenses. 

The asymmetric two-image lenses (those with a large 
magnification ratio between the two bright quasar im- 
ages) are best because for those systems, the mean mag- 
nification of the central image is generally the largest, 
for a given lens galaxy (Mao et al. 2001; Bowman et al. 
2004) . In a highly asymmetric system, the angular sepa- 
ration between the lens galaxy and the unlensed source is 
nearly as large as possible, while still being close enough 
to produce multiple images. This results in a central im- 
age that is as distant as possible from the lens galaxy 
center. More symmetric systems (especially those that 
produce four bright images) have central images located 
closer to the galaxy center, where the surface density is 
larger and the degree of demagnification is consequently 
greater. 

Observing at radio wavelengths is desirable to avoid 
extinction by dust within the interstellar medium of the 
lens galaxy. Attcntuation by plasma effects is also possi- 
ble, but can be minimized by observing at a high enough 
radio frequency (typically > 5 GHz). Moreover, an op- 
tical image could easily be lost in the starlight of the 
lens galaxy, while even a demagnified radio image will 
be brighter than a typical radio-quiet lens galaxy. The 
angular separation between the central image and the 
other images, or the lens galaxy center, is likely to be 
100 mas or less, so at the widely used radio frequencies 
of 1-10 GHz, very-long-baseline interferometry (VLBI) is 
needed. 

The search for central images is a task that is well- 
matched to recent advances in VLBI technologies. New 
high-bandwidth recorders and digital back-ends have in- 
creased typical data rates by a factor of 4 in recent 
years, with even greater improvements expected in the 
next few years. The effect is a considerable sensitivity 
boost for those experiments that can take advantage of 
the increased bandwidth. The European VLBI Network 
now routinely records at 1 Gb s _1 . The U.S. National 
Radio Astronomy Observatory does not yet generally 
deploy recorders with bandwidths this wide, but it re- 
cently began devoting approximately 300 hours per year 
to simultaneous observing with the 10-station Very Long 
Baseline Array (VLBA) and several other large antennas 
such as the phased Very Large Array (VLA), Arecibo, 
the Green Bank Telescope, and the Effelsberg telescope. 
The thermal noise level with this "High Sensitivity Ar- 
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TABLE 1 
Details of the VLBA observations. 



Date 


Image A flux 


density 


Image B 


R.A.a-R.A.b 


dec.A-dcc.B 


Blank field 




point source 


total flux density 






rms 




(mjy) 


(mjy) 


(mjy) 






(/ijy beam -1 ) 


2000 Oct 07 


206.8 


219.2 


9.3 


0'.' 09747 


C'99130 


86 


2000 Oct 30 


204.3 


213.4 


16.2 


0" 09747 


0"99125 


88 


2000 Dec 17 


222.4 


232.3 


12.8 


C'09746 


C'99122 


95 


2000 Dec 18 


202.7 


214.4 


12.3 


0'.' 09747 


/ /99121 


98 


2000 Dec 22 


209.7 


209.7 


13.4 


0'.' 09745 


0'.'99118 


93 


2001 Mar 24 


207.5 


219.5 


18.3 


0'.' 09744 


0'.'99111 


79 



Note. — Details of the individual observing epochs. The beam size varied slightly between epochs, with average values of 3.5 X 1.3 mas, 
and the position angle was always within 1° of zero. All epochs excluded the Brewster antenna due to low elevations. The December 
epochs lacked the North Liberty VLBA antenna, December 17 and December 18 substituted a single Very Large Array (VLA) antenna for 
the Pic Town VLBA antenna. 



ray" (HSA) is often 10-20 times lower than a traditional 
VLBA observation of similar duration. Gravitational 
lens central image searches are perfect projects for sensi- 
tive observations with heterogeneous arrays. The search 
for a faint central image in an otherwise blank region 
is a problem limited mainly by the thermal noise level, 
rather than the need for complete coverage of the Fourier 
plane. The typical radio lens field has two or four bright, 
compact sources (the bright images of the lensed quasar), 
and no extended structure. The bright images provide 
in-beam phase calibration sources. The simplicity of the 
source structure (multiple, isolated point sources) makes 
it easier to calibrate and image the data from the het- 
erogeneous array. 

For all these reasons, the first ELVIS targets are radio- 
loud quasars lensed by a single galaxy to produce two 
bright images with a flux ratio exceeding 5:1. Two lenses 
that fit this description, PMN J1632-0033 (Winn, J. N. 
et al. 2002) and CLASS B1030+074 (Xanthopoulos, E. 
et al. 1998), have already been the subjects of sensitive 
VLBI observations by other groups. In the former case, 
good evidence for a central image was found (Winn et al. 
2004) . We are in the process of obtaining data from long- 
duration VLBA observations, or shorter-duration HSA 
observations, of the other known radio lenses meeting 
these criteria, many of which are taken from CLASS. 

Here we present our first target, PMN J1838-3427. 
This gravitational lens was discovered by Winn, J. N. et 
al. (2000), as part of a survey for radio- loud gravitational 
lenses in the southern sky. It has a two bright images sep- 
arated by l'.'O and with a flux ratio of 14:1. The rcdshift 
of the source quasar is at = 2.78. A spectroscopic red- 
shift for the lens galaxy has not been measured, despite 
several attempts. It is a difficult measurement mainly 
because the optical light of the lens galaxy is blended 
with that of the fainter lens image, and also because the 
system resides in a crowded field at low galactic latitude. 
Winn, J. N. et al. (2000) estimated z h = 0.36 ± 0.08, 
the range of redshifts for which the lens galaxy photom- 
etry is consistent with the fundamental plane relation. 
The source was monitored at 9 GHz with the Australia 
Telescope Compact Array (Winn, J. N. et al. 2004), as 
discussed in Section 3.1. 

3. OBSERVATIONS OF J1838-3427 

We observed PMN J1838-3427 (hereafter, J1838) with 
the ten antennas of the NRAO Very Long Baseline Ar- 



ray (VLBA), on 6 different epochs between 2000 October 
and 2001 March. At each epoch, the duration of the ob- 
servation was 5 hours. We observed right circular polar- 
ization with a central frequency of 8.415 GHz. We used 2 
bit sampling at 16 Msamplcs s _1 for each of 8 channels, 
giving a total data rate of 256 Mb s _1 . Data from anten- 
nas that were observing at an elevation of less than 10 
degrees were excluded (this included all of the data from 
the Brewster antenna). In addition, the North Liberty 
antenna was not in use for two epochs, and was excluded 
from a third epoch for which its data were very noisy 
and degraded the image. Thus, 3 epochs used 9 anten- 
nas and 3 epochs used 8 antennas. In each observation, 
the time on J1838 was 4.5 hours and the estimated ther- 
mal noise limit (based on the collecting area and receiver 
characteristics) was <~ 70 /iJy beam" 1 . 

For each epoch the data were amplitude calibrated in 
AIPS, following standard procedures. Because image A 
of the gravitational lens is fairly bright, with a flux den- 
sity of ~ 200 mjy, we did not need to perform phase 
referencing. The initial phase solution was derived from 
a fringe fit to a point-source model centered at the loca- 
tion of image A, using a solution interval of 2 minutes. 
We then reduced the data with nine self-calibration cy- 
cles, each cycle consisting of imaging, phase-only self- 
calibration with a 0.5-1.5 minute solution interval, imag- 
ing, and phase and amplitude self-calibration with a 15 
minute solution interval. For each new cycle, smaller 
CLEAN components were included. 

The brightest image A appeared as a point source of 
flux ~ 210 mjy, with ~ 10 mjy of extended emission 
to the west, while image B appeared as a point source, 
varying between 9 and 18 mjy. Maps from the first epoch 
are shown in Figure 2, while details of each observing run 
are presented in Table 1. The single epoch 5 GHz VLBA 
map of Winn, J. N. et al. (2000) shows two bright point 
sources at the same positions, with diffuse emission to the 
west of image A. The fraction of the total image A flux 
density in this extended component is ~ 10% at 5 GHz 
and ~ 5% at 8.4 GHz, so the extended emission has a 
steeper spectral index than the point source. 

The flux densities of the two bright images varied be- 
tween the epochs, and the ratio of these flux densities 
also varied. The image A to image B flux density ratio 
in our maps varies from 12 to 24, with an average of 15. 
This ratio averaged 14.6 in VLA observations and 10.6 in 
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Fig. 1. — Radio map of the two images of gravitational lens 
PMN J1838-3427, made with the VLBA at 8 GHz. Data are from 
the first epoch. The restoring beam (40 X 20 mas) was chosen to 
be much larger than the naturally weighted beam in order to show 
both images on a single map. Contours begin at 1.5 mjy beam -1 
and increase by factors of 2. The J2000 radio positions were not 
determined from our observations; they were assumed from earlier 
VLA imaging. The cross marks the location of the lens galaxy 
detected using HST/WFPC2; its position is based on the offset 
from the bright quasar images in the WFPC2 images. The box 
shows the central image search region. 



the the previous 5 GHz VLBA observations (Winn, J. N. 
et al. 2000). The lens images appear to undergo large 
variations in intensity. The two most obvious possible 
explanations are intrinsic source variation, and interstel- 
lar scintillation (see Section 3.1). 

No additional sources of radiation were seen in any 
epoch. We examined a large region between the location 
of the lens galaxy and the bright image B, where a cen- 
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Fig. 2. — (left) Radio map of image A. (right) Radio map of image 
B. Both maps were made from the first epoch of VLBA imaging 
at 8 GHz. Coordinate offsets are from the phase center at (J2000) 
18 h 34 m 28?495, -34°27'41'/60. Contours begin at 300 ^tjy beam" 1 
and increase by factors of 2 (the blank field rms near image B was 
86 fi3j beam -1 ). The synthesized beam of 3.6 X 1.1 mas is shown 
at the lower left. Note the extended emission to the west of the 
main point source in image A. 



I 




♦ 



18 38 28491 28.490 28/489 28.488 28.487 28486 28.485 
RIGHT ASCENSION (J2000) 

Fig. 3. — Radio map of the central image search region, from the 
combined VLBA 8 GHz map of all six epochs. The wedge at right 
shows the grey scale in /tjy beam -1 . The J2000 radio positions 
were not determined from phase referencing; they were assumed 
from earlier VLA imaging. The cross marks the location of the 
lens galaxy detected using HST/WFPC2. Its position is based 
on the offset from the bright quasar images in the WFPC2 images. 
Image B is obvious near the bottom of the map. The central image 
would form at the light travel time maximum located between the 
lens galaxy and the bright image B, so we searched for the central 
image in the boxed area. 



tral image would be expected. 3 We extended this region 
slightly to the other side of the lens galaxy to allow for the 
uncertainty in the optical position. From epoch to epoch 
the noise varied from 79 /iJy beam" 1 to 98 /zJy beam -1 
(Table 1), which is within 20-30% of the expected ther- 
mal noise limit. The small but significant excess over the 
thermal noise limit is expected, given the southerly dec- 
lination of the target. When the six maps were co-added, 

3 Gravitational lens images form at extrema of the time delay 
surface, which is a combination of the geometric delay and the lens 
galaxy potential determined by its density profile. The central im- 
age forms near the lens galaxy potential maximum. The geometric 
delay shifts the overall maximum towards the saddle point image 
B, on the opposite side of the lens galaxy from the source. The 
central image will thus form between the lens galaxy and image B, 
the fainter of the two bright images. 
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Fig. 4. — Histogram of surface brightness values for all pixels 
within the central-image search region. The pixel size is 0.2 mas, 
which is smaller than the typical beam size of 3.5 X 1.3 mas. 
The distribution of surface brightness is approximately Gaussian, 
with a mean of —0.4 ^ijy beam -1 and a standard deviation of 
38 /ijy beam -1 . 



the central-image search region remained blank (see Fig- 
ure 3) and the noise decreased by ~\/6. The distribu- 
tion of surface brightness in the pixels of the final map is 
roughly Gaussian, with a mean of —0.4 pjy beam -1 and 
a standard deviation of 38 /iJy beam" 1 (see Figure 4). 

To quote an upper limit, we considered directly the 
surface brightness distribution and took the 99th per- 
centile to be the 99% limit on the central image flux 
density. This gave an upper limit on the central image 
flux density of 5c < 83 pJy, and a lower limit on the 
magnification ratio of the brightest image to the central 
image of Sa/Sq > 2500, both with 99% confidence. For 
comparison, the 5 GHz VLBA maps of Winn, J. N. et 
al. (2000) reached an rms of 190 /zJy beam -1 , and mea- 
sured a flux density of Sa = 145 mJy for the brightest 
image. Assuming a Gaussian distribution for the blank 
field errors, the old 99% limits were Sc < 470 ^Jy and 
Sa/Sc > 310. 

3.1. Scintillation 

We found the flux densities Sa and Sb of the two bright 
quasar images to vary from epoch to epoch. The ratio 
Sa/Sb also varied significantly, demonstrating that the 
variations are not due only to inconsistencies in the flux 
density scale. Image B showed a higher fractional varia- 
tion than image A (Table 1). This is similar to what was 
observed by Winn, J. N. et al. (2004), who monitored 
this object for 4 months with the Australia Telescope 
Compact Array (ATCA) at 9 GHz. Over the course of 
the campaign, the fractional variation in the flux density 
of image A was 4%, and that of image B was 8%. 

The lens is located at a low galactic latitude (&n = 
— 12?5), and may undergo scintillation due to the Milky 
Way's interstellar medium. The root-mean-squared am- 
plitude due to scintillation is inversely proportional to the 
angular size of the source (Walker 1998, 2001), while the 



angular size of each lens image is proportional to 5 1 / 2 , 
since gravitational lensing conserves surface brightness. 
Thus the rms amplitude of scintillation variations is pro- 
portional to S" 1 / 2 , and the fainter image B would be 
more affected by interstellar scintillation than image A. 
In contrast, intrinsic variability should produce the same 
fractional variation in each component (though the vari- 
ations would appear with time lags due to the geometric 
and Shapiro delays). The greater fractional variation of 
image B that was observed with ATCA, and in our own 
observations, supports the scintillation hypothesis. 

The central image, if it exists, has a much smaller flux 
density and angular size than either of the bright images. 
Taking Sa/Sq > 2500, scintillation would cause frac- 
tional variations in Sc of order unity. No central radio 
source was seen in any individual observation. Scintilla- 
tion, if present, did not magnify the central image above 
the detection limit for any single epoch (~ 200 £tJy). 
In what follows, we use the blank field of the combined 
map to define our upper limit on the central image flux 
density, but we note that scintillation could be a major 
source of systematic error in this determination. 

4. MASS MODELS 

4.1. Cuspy Matter Profiles 

With this new and more stringent upper limit on the 
flux density of the central image, we can restrict the pos- 
sibilities for the central density profile of the J1838 lens 
galaxy. A simple and realistic model for the mass dis- 
tribution of a massive galaxy is a broken power law (see 
Section 1). We adopt the broken power law density pro- 
file of Munoz et al. (2001), in which the surface density 
and the deflection angle are given by analytic expres- 
sions. This profile varies as p oc r~ n at large radii, and 
as p cx r" 7 at small radii, with a break at radius r^: 



> y > r-t (1 + r 2 /V 2 )(™-r>/ 2 ' y > 

Since galactic profiles are approximately isothermal on 
scales of a few kpc (see Section 1) and we have few con- 
straints, we apply the method used by Winn et al. (2003) 
in the analysis of PMN J1632-0033 and fix the outer 
power law to be isothermal (n = 2). We then explore 
the constraints on the break radius r& and inner power 
law index 7. In the limit of a pure isothermal sphere, 
the central image vanishes; thus, as rj, goes to zero or as 
7 goes to 2, the central image flux density approaches 
zero. We use a spherical galaxy model, and account for 
non-sphericity in the profile with an external shear at an 
arbitrary position angle. The 10 free parameters of our 
model are the positions of the lens galaxy and the source, 
(xq, VGi xs, Us), the source flux (5s), the mass parameter 
(b = (27rrb/9o)/(^cr)) 4 the shear and its position angle 
(s,8 s ), the break radius (r b ) and the inner power law 
index (7). 

The data provide 9 observables: the two sky coordi- 
nates for each of the lens galaxy, bright image A and 
bright image B; the flux densities for image A and image 

4 The quantity E cr = (c 2 j AttG){D s / (D L D LS )) is the critical 
surface density for lensing. The quantities D^, Dg and D^g arc 
the angular diameter distances from the observer to the lens, from 
the observer to the source and from the lens to the source, respec- 
tively. 
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TABLE 2 

Lens model constraints. 

Observable Value 
R.A.a O'.'OOO ± 0'.'003 

decl. A O'.'OOO ± 0'.'003 

R.A.b -0"9912 ± 0'.'003 
dccl.B -0"0975 ± 0'.'003 
R.A.g -0'.'911±0'.'006 
decl. G -0'.'085 ± 0'.'006 
S A 208.3 ± 41.8 mjy 

S B 13.7 ± 2.7 mjy 

S c 0-000 ± 0-038 mjy 
Note. — Constraints on the lens models, with 1 cr errors. All 
positions arc offsets from image A. Galaxy postions are from the 
WFPC2 photometry in Winn, J. N. et al. (2000), other constraints 
are from the VLBA data. 



B (Sa, <5b); and an upper limit on the flux density of the 
central image C (Sc)- The galaxy coordinates and their 
uncertainties were taken from the Hubble Space Tele- 
scope Wide Field Planetary Camera 2 (WFPC2) pho- 
tometry by Winn, J. N. et al. (2000). The radio coordi- 
nates and flux densities for images A and B were mea- 
sured in each of the six VLBA maps and the mean val- 
ues of these observables were taken as model constraints. 
We considered only the bright peak of image A, and ig- 
nored the diffuse emission to the west, as we consider 
only the point source magnifications. Since both the 
positions and fluxes of the bright images may be per- 
turbed by substructure in the lens galaxy, we adopted 
errors of 3 mas in the positions and 20% in the flux den- 
sities (cf. Mao & Schneider (1998); Keeton & Kochanck 
(1998); Kochanek et al. (2004)). The constraint on S c 
was taken from the rms noise in a large region between 
the lens galaxy position and image B in the combined 
VLBA image. The distribution of surface brightness in 
the pixels of this blank region was approximately Gaus- 
sian, with a mean of —0.4 LtJy beam -1 and a standard 
deviation of 38 /iJy beam" 1 (see Figure 4). Wc took 
this to be a null detection of the central image with a 
1 a error of 38 /zJy beam, and set Sc = ± 38 /iJy 
for the modelling. These constraints are summarized in 
Table 2. All of the observables were assumed to obey 
Gaussian statistics, and 1 a error bars are quoted. 

We model the lens with the gravlens package by Keeton 
(2001). For the conversion from angular units to phys- 
ical units, we assume zl = 0.36 and a ACDM cosmol- 
ogy with H = 70 km s" 1 Mpc" 1 ,^ = 0.3, fi A = 0.7. 
With 10 free parameters and 9 measurements, one would 
expect a one-dimensional locus of allowed points in pa- 
rameter space. In fact, because the observable Sc is an 
upper limit, the data define a one-dimensional boundary 
between an allowed region and a disallowed region in pa- 
rameter space. In Figure 5, we plot the projection of this 
plane in the n,, 7 parameter space that characterizes the 
inner density profile. Values of the other parameters are 
given in Table 3 for the isothermal case (rt, = 2 or 7 = 0), 
and these values change by < 10% over the region of r^, 7 
parameter space we investigated. 

For a sufficiently small break radius r^ or sufficiently 
steep inner index 7, the model fits the data perfectly. As 
rb increases or 7 decreases, the fits become progressively 
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Fig. 5.— Models for the lens galaxy of PMN J1838-3427, using 
the new data presented in this paper (solid line) and the previous 
data from Winn, J. N. et al. (2000) (dashed line). The matter 
profile is isothermal (p oc r~ 2 ) at large radius, and goes as p oc r ' 
at small radius, with the break at radius r^. The inner power law 
index and break radius are allowed to vary. The lines mark the 
\ 2 = 5.99 boundary between the allowed region (below and to the 
right) and the disallowed region (above and to the left). The break 
radius in parsecs is calculated for zl = 0.36. 



TABLE 3 

Lens model parameters for the isothermal case. 



Parameter 


Value 


b = E 


(X'528 


XQ 


-0'.'0850 


VG 


-0'.'9109 


x s 


-0'.'0187 


ys 


-(X'4747 


s s 


80.3 mjy 


s 


0.0649 


8s 


-18?1 



Note. — The best fit values of the lens model parameters for 
the isothermal case (when r-j, = or 7 = 2). 



worse and \ 2 becomes large. Most of the disagreement is 
due to an unacceptably large Sc- Near the boundary of 
the allowed region the best fits give slight variations in 
the calculated values for the other observables, but these 
are generally less than the 1 a errors. The 95% confidence 
region (\ 2 < 5.99 for two free parameters, cf. Bevington 
& Robinson (2003)) lies below and to the right of the 
solid line. The matter distribution is nearly isothermal: 
a flat core with 7 = must be very small (rt < 1.5mas w 
5 pc), and a large core with rt = lOOmas s=s 500 pc has a 
steep inner profile (7 > 1.92). 

For comparison, Figure 5 includes the region of pa- 
rameter space allowed by the less stringent observational 
constraints of Winn, J. N. et al. (2000). The new data 
reduce the maximum size of a flat core by a factor of 
^2.5, and increase the minimum inner index of a large 
core by ~ 0.05. 

4.2. Cuspy Matter Profiles with Black Holes 

Most elliptical galaxies host super-massive black holes 
at their centers (Kormendy & Richstone 1995; Magor- 
rian, J. et al. 1998): this additional point mass steepens 
the overall central profile and may demagnify the central 
image (Mao et al. 2001). A realistic galaxy model should 
include a central black hole. As seen in the previous 
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section, even an interesting smooth profile is undercon- 
strained. To avoid adding extra degrees of freedom to the 
model, we add a black hole with a fixed mass given by 
the M — a relation (Ferrarese & Merritt 2000; Gebhardt, 
K. ct al. 2000; Tremaine, S. et al. 2002), and determine 
the resulting effect. 

For an early-type lens galaxy, most of the mass will be 
in the bulge component and so the total mass responsible 
for the lensing is a good estimate of the overall velocity 
dispersion. To estimate ct, we determine the best-fitting 
value of the Einstein radius using a model consisting of 
an isothermal sphere and external shear, and then use 
the relation f?E = 47t(£)ls / 'Ds)(a 2 / ' <?) ■ For the broken 
power law model, the mass normalization parameter b — 
(2TrrbPo)/(Ecr) is equivalent to #e and this parameter 
varies by < 10% over the allowed region in the broken 
power law models, so the isothermal #e is a good estimate 
of the true velocity dispersion even if the true density 
profile is not exactly isothermal. 

Given this estimate for ct, we calculate the black hole 
mass Mbh using the relation from Tremaine, S. et al. 
(2002), 



M bh = 1.35 x 10 8 M ( 



200 km s 



-) 



4.02 



(2) 



For J1838, E = 0'.'523, ct = 154 km s" 1 and M bh = 
4.7 x 10 7 M Q . The scatter in the M — a relation is 
0.3 dex (Tremaine, S. et al. 2002), so a plausible range 
for the black hole mass is 2.5 x 10 7 - 1.0 x 10 8 M . 

We rerun the calculations from Section 4.1, but this 
time include a point mass fixed at the center of the lens 
galaxy. We consider a black hole with mass 5 x 10 7 M , 
in the center of the expected range, and a black hole with 
the mass 10 8 M , at the top of the expected range. The 
black hole either eliminates the central image completely, 
leaving only the two bright images, or creates two faint 
central images, the core and SMBH images, for a total 
of four images. The allowed regions, again defined by 
X 2 < 5.99, are shown in Figure 6. Below and to the right 
of the allowed lines there are no central images, or both 
the core and SMBH images are sufficiently faint that they 
would not have been detected in our observations. The 
long-dashed line is for the 5 x 10 7 M black hole, the 
short-dashed line is for the 10 8 M black hole. 

The separation of the core and SMBH images is 2 — 
5 mas, so these might appear as a single unresolved com- 
ponent if aligned with the long axis of the beam. The 
core image flux density is greater than the SMBH image 
flux density by a factor of 15 — 100, so the core image 
flux density makes a much larger contribution to x 2 ■ F° r 
models near the boundary of the allowed region, com- 
bining both central images into a single unresolved com- 
ponent would raise \ 2 by < 0.6, and the allowed region 
would be only slightly more constrained. 

The SMBH allows the galaxy profile to be shallower, to 
some extent. For fairly flat cores (7 < 1.5) the maximum 
break radius increases by about 50% for the 5 x 10 7 M 
black hole, and by a factor of 2.5 — 3 for the 10 8 Mq black 
hole. For larger, cuspy cores, the minimum inner slope 
is reduced by <~ 0.1 for the 5 x 10 7 M Q black hole and 
by ~ 0.2 for the 10 s M black hole. While the central 
smooth profile remains close to isothermal, it is allowed 
to be somewhat flatter. The cusp in the profile from 
the black hole makes the central mass profile sufficiently 
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Fig. 6.— Models for the lens galaxy of PMN J1838-3427, us- 
ing the new data presented in this paper and including the effects 
of a super-massive black hole. The matter profile is isothermal 
(p oc r~ 2 ) at large radius, and goes as p <x r ' at small radius, 
with the break at radius r;,. To the smooth matter profile we add 
no black hole (solid line), a 5 X 10 7 Mq black hole (long-dashed 
line), or a 10 8 Mq black hole (short-dashed line). These black hole 
masses are the central and largest values expected, given the veloc- 
ity dispersion of the galaxy. The inner power law index and break 
radius arc allowed to vary. The lines mark the \ 2 = 5.99 boundary 
between the allowed region (below and to the right) and the disal- 
lowed region (above and to the left). The break radius in parsccs 
is calculated for zl = 0.36. When no black hole is present, a cen- 
tral image always forms within the parameter ranges plotted. In 
the allowed region, the central image is below our detection limit. 
When a black hole is present, the central image is cither eliminated 
or split into a core image and SMBH image. In the allowed region, 
either there are no central images or both central images have flux 
densities below our detection limit. 



steep even if the smooth component is shallower. 

4.3. Comparison with Other Central-Image Modelling 

Winn et al. (2003) modelled PMN J1632-0033, the lens 
with a confirmed central image, using the same cuspy 
profile as this work. Given that the central image is 
detected with a definite flux, the allowed region in the 
break radius and inner index space becomes a narrow 
band (see Winn et al. (2003), Figure 6). A flat core has a 
fixed small size, and a larger core is forced to a particular 
value of the inner index. The allowed band for J 1632 is 
similar to the boundary of the allowed region for J 1838: 
a central image of definite flux density selects a narrow 
band giving that value for Sc, while an upper limit on 
the central flux density selects a half plane bounded by 
the band for the maximum allowed Sc- 

Other central image modelling has focused on single 
power law models. Winn ct al. (2003) also considered 
a single power law model p cx r~ 7 for the density pro- 
file of J1632, finding 7 = 1.91 ± 0.02. This modelling 
did not include a central black hole. Rusin & Ma (2001) 
modelled the gravitational lenses CLASS B0739+366 and 
B1030+074 using a single power law model for the pro- 
jected surface density S(i?) oc R~P, both with and with- 
out a central black hole with a mass inferred from the 
M — a relation. Their only constraint was the limit on 
the bright image to central image magnification ratio, the 
observable which dominates \ 2 m our modelling. Their 
limits on the smooth profile were > 0.85 for B0739+366 
and (3 > 0.91 for B1030+074, while the respective lim- 
its were (3 > 0.84 and (3 > 0.83 when a black hole was 
included. 
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Modelling J1838 as a single power law in the den- 
sity p oc r~ 7 , we find 7 > 1.93 for the smooth pro- 
file, 7 > 1.86 when we include a central black hole of 
5 x 10 7 M Q and 7 > 1.82 when we include a central 
black hole of 10 8 M© (the central and largest values ex- 
pected from the velocity dispersion). Again, we note that 
a surface density index (3 is equivalent to a density index 
7 = (3 + 1, so these 4 lenses are similar. The profile is 
slightly shallower than isothermal for J1632-0033, and 
the profiles are slightly shallower than or steeper than 
isothermal for J1838-3427, B0739+366 and B1030+074. 

5. CONCLUSIONS 

We have presented new VLBA data for the gravita- 
tional lens PMN J1838-3427, and set a stringent upper 
limit of 83 /xJy (99% confidence) on the observed flux 
density of the central image. The corresponding limit 
on the brightest image to central image magnification 
ratio is Sa/Sq > 2500. This improves on the limits ob- 
served by Winn, J. N. et al. (2000) by a factor of 8. The 
bright images have different flux densities and different 
flux density ratios in our six epochs, confirming the fluc- 
tuations seen by Winn et al. (2003). The larger fractional 
variation in image B relative to image A is suggestive of 
scintillation. The possibility that the central image also 
undergoes scintillation is a source of systematic uncer- 
tainty in the interpretation of the upper limit on Sq- 

Lens galaxies at z < 0.5 have isothermal profiles on 
scales of a few kiloparsecs (Rusin et al. 2003; Koopmans 
et al. 2006), and we find that the lens galaxy in J1838 
maintains this steep profile to within a few tens of par- 
sees of its center. The lenses J1632-0033, B0739+366 
and B1030+074 also have central density profiles close 
to isothermal (Winn et al. 2003; Rusin & Ma 2001). 
These four lens galaxies resemble the nearby cusp galax- 
ies in the photometric observations of Lauer, T. R. et al. 



(1995). Galaxy simulations which track both dark mat- 
ter and baryons predict isothermal profiles at galactic 
radii > 1 kpc (see Section 1), and central image con- 
straints from these galaxies show these profiles contin- 
uing to smaller radii. A black hole with a mass esti- 
mated from the M — a relation weakens the constraints 
somewhat, allowing the central profile to be a little more 
shallow. 

As we add more galaxies to the ELVIS sample, we will 
have a larger sample of gravitational lenses with either 
central images or very good upper limits on the central 
image flux density. Central images are very sensitive 
to the central matter distributions that are just shal- 
lower than the isothermal profiles of many local galaxies 
and of distant galaxies at larger radii. With more well- 
constrained central matter profiles, we can determine if 
most early type galaxies at z ~ 0.5 have cuspy central 
profiles, or if some have shallow cores. As well as the 
central image from the galaxy core, we may detect cen- 
tral image pairs and directly measure black hole masses 
in distant galaxies. 
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